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Shuttle  Contamination  Modeling: 
Evolution  of  Ionized  Shuttle  Exhaust 


1.  INTRODUCTION 

One  important  aspect  of  the  interaction  of  Shuttle  with  its  environment  is  the  evolu¬ 
tion  of  contaminants  and  exhaust  of  Shuttle  origin.  The  purpose  of  this  paper  is  to  provide 
an  order  of  magnitude  survey  of  the  interaction  of  the  exhaust  with  the  ambient  environ¬ 
ment.  The  topics  considered  include  the  collisional  and  plasma  behavior  of  the  exhaust 
and  the  interaction  of  ionized  exhaust  with  the  ambient  environment  and  electromagnetic 
fields.  While  the  sources  of  ionization  are  not  treated,  the  subsequent  evolution  is  treated 
in  a  certain  amount  of  detail. 

The  focus  of  the  paper  is  the  evolution  of  ionized  exhaust.  In  particular,  the  goal  is  to 
specify  the  trajectory  of  the  ionized  exhaust  and  to  determine  the  magnitude  of  the  effects 
of  ambient  electric  and  magnetic  fields.  We  will  examine  the  behavior  of  exhausts  with  low 
and  high  fractional  ionization  and  discuss  the  qualitative  differences  in  their  evolution. 

The  fundamental  conclusions  are  that:  (1)  Low  density  ions  expand  with  the  neutral  ex¬ 
haust  for  about  100  m  where  they  go  through  a  rapid  transition  to  collisionlcss  behavior 
and  then  arc  swept  away  (in  the  Shuttle  frame  of  reference)  by  the  geomagnetic  field,  and 
(2)  High  density  ions  are  not  affected  by  the  magnetic  field  and  arc  expected  to  expand 
radially  away  from  the  Shuttle;  neither  the  ambient  atmosphere  nor  the  geomagnetic  field 
can  significantly  change  this  behavior.  The  details  of  the  evolution  can  be  affected  by  low 
level  (  =  0.1  mV/m)  parallel  electric  fields. 

The  motive  of  the  paper  is  to  identify  some  of  the  physical  issues  that  must  be  dealt 
with  in  an  effort  to  model  contamination  processes  on  Shuttle  and  other  large  spacecraft. 
The  release  of  Shuttle  exhaust  has  several  possible  ways  of  contributing  to  the  contami¬ 
nant  environment.  Among  them  are  (1)  direct  contamination  of  Shuttle  surfaces  and 
sensors;  (2)  obscuration  of  the  field  of  view  of  sensors,  principally  by  infrared  absorption 
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and  (3)  optical  contamination  by  glow  associated  with  excited  or  ionized  atoms.  While  our 
discussion  could  be  applied  to  all  three  categories,  we  are  directly  concerned  only  with  the 
third.  The  nature  of  our  inquiry  is:  What  is  the  spatial  evolution  of  Shuttle  exhaust,  as¬ 
suming  that  it  is  partially  ionized?  The  intent  is  to  be  able  to  outline  the  physical  parame¬ 
ters  necessary  to  determine  the  evolution  and  motion  of  an  ion  cloud. 

Despite  the  fact  that  we  are  able  to  analyze  the  motion  to  a  certain  extent,  the  impor¬ 
tance  of  the  results  is  that  they  point  out  the  physical  basis  for  modeling  the  behavior  and 
evolution  of  the  exhaust.  As  such,  the  present  paper  should  be  viewed  not  as  a  complete 
analysis  of  the  problem,  which,  in  fact,  it  is  not  intended  to  be,  but  as  a  starting  point  for 
such  an  analysis.  We  believe  that  the  proper  treatment  of  such  phenomena  as  Shuttle 
exhaust  and  contamination  requires  computer  simulation.  Yet  much  of  the  behavior  is 
quantifiable  using  analytical  methods.  It  is  our  feeling  that  aspects  of  the  behavior  should 
be  incorporated  into  the  computer  program  where  possible,  rather  than  modeled  by  purely 
numerical  methods. 

The  text  is  organized  into  separate  sections.  The  plan  of  each  section  is  to  give  (l)a 
sequence  of  order-of-magnitude  estimates  designed  to  focus  on  one  aspect  of  the  problem, 
followed  by  (2)  a  discussion  designed  to  help  develop  a  coherent  argument.  A  separate 
discussion  section  at  the  end  of  the  paper  draws  from  the  conclusions  of  each  section. 

2.  TIIK  KXHAI  SI 

The  purpose  of  this  section  is  to  address  the  neutral  properties  of  the  Shuttle  ex¬ 
haust:  composition,  fluid  properties,  and  eollisional  behavior  of  both  neutrals  and  ions.  In 
addition,  we  discuss  simple  properties  of  the  eollisional  interaction  of  the  exhaust  with  the 
neutral  atmosphere.  We  will  be  primarily  interested  in  the  distances  from  the  Shuttle  at 
which  transitions  from  one  type  of  behavior  to  another  occur  (for  example,  from  eollisional 
to  collisionless  behavior). 

The  molecular  composition  of  the  Shuttle  RCS  exhaust  (Bareiss  et  al,‘  p.  E-15)  is 
given  in  Table  1: 


Table  1.  Chemical  Properties  of  Shuttle  RCS  Engine  Exhaust 


Ionization 

Molecule 

Atomic  Number 

Mole  Fraction 

Energy(eV) 

y  =  !■!' 

IDO 

18 

0.33 

12.6 

1.3 

N, 

28 

0.31 

15.58 

1.4 

CO, 

44 

0.036 

13.8 

1.3 

II, 

2 

0.17 

1.4 

CO 

28 

0.13 

14.0 

1.4 

1.  Bareiss,  L.E.,  Jarossy,  F.J.,  Pizzicaroli,  J.C.,  and  Owen,  N.L.  (198D  Shuttle/Payload 
Contamination  Evaluation  Program:  The  SPACE  Computer  Program  User’s  Manual , 
Martin  Marietta  Corp.,  Denver. 
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From  Table  1,  it  follows  that  the  mean  molecular  weight  is  A  =  16.  For  purposes  of  esti¬ 
mating  orders  of  magnitude,  we  will  use  A  =  16.  y  is  the  ratio  of  specific  heats  for  each 
molecular  species.  The  mean  value  of  y  is  slightly  less  than  1.4;  we  will  use  1.4  as  a  typical 
value. 

The  exhaust  of  the  thrusters  has  a  velocity  V  =*  3  km/scc  relative  to  the  orbiter  and  a 
temperature  T  =  3000° K;  the  flux  of  the  exhaust  is  given  by 


|)V  =  -  ..  -cos  (0.01260)^  gni/(cm^sec)  for  0  <  64°  (1) 

r- 


,»V  =  -fexp  -  0.035(0 -64)  gm/(cm^sec)  for  0  >  64° 


where  0  is  in  degrees  and  r  is  in  cm.2  The  flux  is  plotted  in  Figure  1. 

To  obtain  order-of-magnitude  estimates  for  the  fluid  and  collisional  properties  of  the 
exhaust,  we  will  assume  that  the  exhaust  expands  at  constant  velocity: 


V  =  3  km/sec  =  constant. 


From  the  mass  flux  of  the  neutral  exhaust 


(2) 


,,V  =  gm/lcm^scc) 


one  can  then  compute  the  mass  density 


„  _  1350  _  4.5  x  10- 3  m3 
l»  -  Vr2  = - r2 - gm/cm° 


and  the  number  density 


(3) 


(4) 
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Figure  2.  Characteristic  Behavior  of  Mach  Number,  Bulk  Velocity,  Sound  Speed, 
Temperature,  and  Number  density  as  Functions  of  Distance  From  Nozzle 


More  detailed  results  are  shown  in  Figure  2.  These  results,  for  y  =  1.4,  were  computed 
for  adiabatic  expansion  through  a  dc  Laval  nozzle  of  4-cm  diameter.  Because  the  de  Laval 
model  used  is  essentially  one-dimensional,  the  results  are  only  qualitative,  but  they  do 
indicate  the  characteristic  variations  associated  with  adiabatic  expansion.  The  variables 
plotted  are  the  velocity  V,  the  sound  speed  cs  =  \  ykT/m,  the  Mach  number  V  cs,  the 
temperature  T,  and  the  number  density  n.  They  are  all  normalized  to  their  values  at  the 
supersonic  transition  at  the  throat  of  the  nozzle.  They  are  plotted  as  functions  of  distance 
from  the  throat. 


it  =  4rrx  10  ^cm- 


(8) 


is 


X  =  3.0x10  ‘Ar~cm. 


(9) 


The  distance  at  which  the  mean  free  path  is  equal  to  the  scale  height  L  ~  £  i-s 


r  = 


(10) 


For  A  =  16,  this  is  about  1  km.  This  suggests  that  the  strongest  interaction  between 
ambient  molecules  (including  atoms  and  ions)  and  exhaust  molecules  occurs  at  about 
1  km.  Outside  this  distance,  the  exhaust  is  so  tenuous  that  collisions  with  ambients  are 
rare.  Inside  this  distance,  the  collisional  mean  free  path  is  small  enough  that  ambients  can 
penetrate  only  with  difficulty.  This  is  clearly  only  an  order-of-magnitude  estimate;  the 
details  of  the  penetration  depend  on  the  relative  velocities  and  the  angle  at  which  the 
ambients  enter  the  exhaust.  In  particular,  ambient  flux  impinging  on  the  sides  of  the 
exhaust  cone  will  have  a  strong  interaction  over  about  a  mean  free  path,  which  is  short 
near  the  nozzle.  Similar  arguments  can  be  developed  for  the  exhaust  outside  the  main 
plume.  There  the  mass  flux  is 


,,V  ~ 


35 


r- 


gm 

(cm'-’sec) 


(Hi 


and  the  distance  where  the  mean  free  path  is  equal  to  the  scale  height  is  about  30m. 

An  important  consequence  of  the  collisional  behavior  of  the  exhaust  is  that  ions  in  the 
exhaust  are  necessarily  swept  out  with  the  exhaust  until  a  transition  to  a  collisionless 
regime  is  reached.  To  anticipate  the  results  of  the  next  section,  the  gyroperiod  of  an  ion 
with  a  velocity  equal  to  that  of  the  thruster  is  about  1.6  x  10  A  sec  (for  singly  ionized 
ions).  The  mean  free  time  of  a  thermal  ion  is 


1 

(nowi 


(12) 


where  w  is  the  thermal  speed  of  the  exhaust: 


2kT 

m 


(13) 


and  k  is  Boltzmann’s  constant.  The  distance  at  which  the  mean  free  time  is  equal  to  the 
gyroperiod  is 


r 


400 
\  M 


(14i 


where  M  =  V,cs  —  V/w  is  the  Mach  number.  Up  to  this  distance,  the  ion  suffers  many 
collisions  per  gyroperiod  and  is  swept  out  by  the  exhaust.  The  implication  of  multiple 
collisions  is  that  the  ions  are  thermalized  in  the  frame  of  reference  of  the  exhaust;  in  the 
Shuttle  frame,  an  ion  moves  with  the  bulk  speed  of  the  exhaust  plus  its  thermal  motion. 
After  this  point,  the  ion  suffers  few  collisions  per  gyroradius,  and  ions  in  the  exhaust 
follow  smooth  trajectories  in  the  local  electric  and  magnetic  fields.  The  subsequent  behav¬ 
ior  of  the  ions  depends  on  the  properties  of  the  plasma  rather  than  the  collisional  behavior 
of  the  neutral  gas. 

The  transition  to  collisionless  behavior  is  rather  abrupt  if  the  exhaust  has  cooled 
rapidly.  Repeating  the  above  argument  while  using  the  fact  that  the  mean  free  path  in  the 
Shuttle  frame  depends  on  velocities  ranging  from  roughly  V-w  to  V  + w,  one  sees  that  the 
range  of  distances  implied  by  Eq.  (14)  is 


r 


400 

\  (M  ±  1 1 


m 


so  that  the  thickness  of  the  transition  region  for  large  M  is  on  the  order  of 


Ar  =  400M  3“  m 


(15) 


(10) 


From  Figure  2,  the  Mach  number  (on  the  assumption  of  adiabatic  expansion)  beyond  a 
few  meters  is  about  10  for  y  =  1.4,  so  the  transition  distance  is  about  100m.  The  thick¬ 
ness  of  the  transition  region  on  this  estimate  is  therefore  of  the  order  of  100/M  —  10  m. 
Since  any  observed  transition  is  unlikely  to  be  less  than  a  gyroradius  thick,  a  better  esti¬ 
mate  may  be  Am  for  heavy'  (singly  ionized)  ions.  In  any  event,  it  should  be  noted  that  the 
maximum  extent  of  the  region  occupied  by  ions  is  the  transition  distance  plus  about  one 
gyroradius  (for  a  total  -  100  to  110m)  if  the  ions  do  not  continue  to  move  along  with  the 
neutral  exhaust. 

The  Shuttle  RCS  exhaust  is  delivered  in  bursts  of  80-300  m  sec  and  the  total  mass  of 
each  burst  is  —  142gm.  The  number  of  particles  in  each  burst  is  (for  A  =  10) 


For  purposes  of  comparison,  the  number  density  of  the  neutral  ionosphere  is 


n  -  10®  cm  ^ 


(18) 


so  the  mass  density  is 


i>  -  3x  10  ^gm.'cm^ 


(19) 


The  total  mass  in  lknv^  of  ambient  ionosphere  is  about  3gm.  As  a  result,  the  exhaust 
dominates  the  neutral  ionosphere  near  the  Shuttle.  The  distance  at  which  the  neutral 
ionosphere  begins  to  dominate  is  where  the  exhaust  mass  density  has  dropped  to  the 
neutral  ionospheric  mass  density: 


r  —  4  km 


(20) 


The  ionization  state  of  the  exhaust  is  not  well  known.  Results  of  the  computer  pro¬ 
gram  CONTAM  III,  which  models  the  combustion  process  inside  the  rocket  chamber  but 
docs  not  account  for  recombination  processes  during  exit,  show  the  number  of  neutral 
particles  ejected  to  be  -  10^5  and  the  number  of  ejected  ions  to  be  -  lO1^  for  a  fractional 
ionization  of  -  10~8.:l  This  result  is  consistent  with  estimates  based  on  the  Saha  equa¬ 
tion  (sec  below,  this  section).  Cooling  of  the  exhaust  and  recombination  processes  may 
decrease  the  ionization  as  the  exhaust  expands.  However,  for  reference  purposes,  let  us 
take  this  figure  as  representative  of  the  state  of  ibe  exhaust  outside  the  orbiter.  The  total 
number  of  ions  released  per  burst  is 


Nj  =  1016.  (21) 

A  crude  estimate  of  the  minimum  volume  occupied  by  the  ions  in  one  burst  for  particles 
injected  along  the  field  is 


E  =  Trp^VAt 


(22) 


where  p  =  10m  is  the  ion  gyroradius  (based  on  the  ion  thermal  speed)  and  At  is  the  time 


3.  Pickett,  Jolene  S.,  Murphy,  Gerald  B.,  Kurth,  William  S.,  and  Goertz,  Christopher  K. 
(1985)  Effects  of  chemical  releases  by  the  STS  3  orbiter  on  the  ionosphere, 

J.  Geophys.  Res.  90:3487-3497. 


duration  of  the  burst.  This  assumes  that  the  burst  occupies  a  column  of  length  VAt  and 
cross-section  irp^.  The  corresponding  maximum  ion  density  is 


n-  =  _i?L_ 

‘  T7p*VAt 

=  105cm'3. 


(23) 


If  the  exhaust  is  injected  across  the  magnetic  field,  the  maximum  number  density  is 


= _ Nj _ 

1  2  x  irp^wAt 


(24) 


where  p  =  20  m  is  the  gyroradius  based  on  the  velocity,  and  the  extra  factor  of  two 
accounts  for  the  spreading  of  the  beam  in  both  directions  along  the  field,  so  that 


nj  —  10“*  cm  3. 


(25) 


These  estimates  show  that  the  maximum  ion  number  density  from  the  exhaust  is  of  the 
same  order  of  magnitude  as  the  ambient  ion  density  (nj  -  104cm-3). 

While  these  estimates  establish  maximum  number  densities,  they  are  probably  too 
high  for  two  reasons:  (1)  The  temperature  is  poorly  known,  and  our  estimate  may  be  too 
high  because  of  adiabatic  cooling;  (2)  The  ions  will  be  spread  over  a  large  surface  at  the 
transition  to  collisionless  behavior.  The  effect  of  lowering  the  temperature  is  to  lower  the 
number  of  ions  in  the  exhaust  by  up  to  many  orders  of  magnitude  because  of  the  tempera 
ture  dependence  in  the  Saha  equation.  The  effect  of  spreading  the  ions  over  the  transition 
layer  is  to  lower  the  ion  density  by  about  one  order  of  magnitude. 

The  largest  uncertainty  in  establishing  the  intrinsic  ionization  state  of  the  exhaust  is 
the  uncertainty  in  the  temperature  of  the  exhaust  as  as  function  of  distance.  Even  moder¬ 
ate  uncertainties  in  the  temperature  cause  large  uncertainties  in  the  fractional  ionization. 
This  sensitivity  follows  from  the  temperature  dependence  of  the  Saha  equation 

r«  •5O0'®(S)«p(^j  .20, 


where  the  fractional  ionization  is  defined  as  the  ratio  of  the  number  density  of  ions  to  that 
of  neutrals,  nn  is  the  number  of  neutrals  per  cubic  meter,  and  I  is  the  ionization  energy 
(ef.  Chen,'  p.  1).  The  striking  feature  of  Eq.  (26)  is  that  if  the  temperature  is  constant, 


4.  Chen,  F.  (1984  Introduction  to  Plasma  Physics,  Plenum,  New  York. 


the  fractional  ionization  rises  with  distance  because  of  the  inverse  dependence  on  the 
neutral  density  while  if  the  temperature  decreases  adiabatically,  the  fractional  ionization 
is  essentially  zero  past  a  few  meters  from  the  aperture. 

This  extreme  sensitivity  to  the  details  of  the  temperature  profile  leads  to  the  following 
observation:  The  nominal  composition  of  the  exhaust  given  in  Table  1  may  be  useless  to 
determine  the  ionization  state  of  the  exhaust.  Low  level  (for  example,  1  part  in  10®  or 
108>  contaminants  with  low  ionization  potentials,  such  as  metals,  would  lead  to  much 
higher  ionization  levels  than  the  nominal  constituents,  assuming  that  the  exhaust  expands 
adiabatically. 

A  number  of  important  results  follow  from  these  order-of-magnitude  estimates: 

(1)  The  ions  in  the  exhaust  will  be  swept  out  with  the  neutral  exhaust  for  about 

100  m.  At  about  100m,  there  is  a  transition  to  collisionless  behavior;  the  transition  region 
is  on  the  order  of  10  m  thick. 

(2)  The  strongest  collisional  interaction  with  the  ambient  atmosphere  is  about  1  km 
from  the  Shuttle.  This  is  the  region  where  neutral  collisions  would  be  most  likely  to  cause 
ionization,  if  ionization  is  energetically  possible. 

<3i  The  ambient  ionosphere  cannot  have  much  efTect  on  the  trajectory  of  the  neutral 
exhaust.  It  simply  does  not  have  enough  momentum  to  stop  more  than  an  extremely  small 
fraction  of  the  exhaust.  As  a  result,  the  neutral  part  of  the  exhaust  can  be  expected  to 
travel  unimpeded  through  the  ambient  environment  for  a  distance  on  the  order  of  4  km. 
Even  at  this  distance,  the  neutral  exhaust  cannot  be  deflected  by  the  ambient  atmosphere 
because  the  mean  free  path  is  too  large.  This  behavior  could  be  changed  by  such  processes 
as  charge  exchange,  which  are  not  considered  here. 

(4 1  Reasonably  high  ionization  states  are  obtained  if  the  exhaust  expands  isother- 
mally;  however,  the  ionization  state  is  so  temperature-dependent  that  adiabatic  expansion 
means  that  the  nominal  exhaust  molecules  are  not  ionized  past  a  few  meters  from  the 
aperture. 

.».  AMIIIENT  ELKCTiUC  AM)  MAGNETIC  FIELDS 

An  ionized  exhaust  evolves  under  the  influence  of  ambient  electric  and  magnetic 
fields.  The  magnitude  of  the  effects  depends  not  only  on  the  strength  of  the  ambient  fields, 
but  the  extent  to  which  the  plasma  can  shield  itself  from  them.  The  purpose  of  this  sec¬ 
tion  is  to  determine  the  order  of  magnitude  of  the  effects  of  the  ambient  fields  on  the 
evolution  of  Shuttle  exhaust.  The  fundamental  assumption  is  that  the  plasma  responds  to 
the  fields  as  a  collection  of  single  particles.  As  a  result,  we  ignore  shielding  and  polariza¬ 
tion;  they  will  be  discussed  in  a  separate  section. 

Typical  Shuttle  altitudes  are  300  km.  Let  us  take  a  typical  geomagnetic  L-shell  to  be 
L  =  1.2  and  the  magnetic  field  strength  to  be  B  =  0.4  G. 

For  this  magnetic  field,  the  gyrofrequency  of  ions  is 

o  -  ZeB 

1  (AmMc)  12,1 

=  3800  (?)  sec-1 
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where  Z  represents  the  ionic  charge  state,  e  is  the  proton  charge,  and  mu  is  the  mass  of 
hydrogen.  For  A  =  16  and  Z  =  1, 


11]  =  240  sec  1 


(28) 


The  ion  gyroradius  for  the  propulsion  velocity  V,  that  is,  for  the  exhaust  directed 
across  the  field,  is 


AmncV 

i>l  ~  TZcBr 


(29) 


or,  for  A  =  16  and  Z  =  1, 


l>]  —  12  m. 


This  estimate  is  appropriate  for  ions  which  have  been  driven  by  collisions  to  a  collisionless 
region.  The  electron  gyroradius  is 


mt.cv 

'’e  “  eB 


V  eB 


(31) 


=  8.4\  t  (eV)  cm 


where  i  is  the  electron  energy  expressed  in  electron-volts.  The  gyroradius  for  electrons 
with  energies  up  to  lOOeV  is  only  about  10cm.  For  the  gyroradius  to  be  as  large  as  10  m, 
the  electron  energy  would  have  to  be  lOkcV. 

The  large-scale  behavior  of  charged  particles  is  determined  by  electric  and  magnetic 
forces  through  the  Lorentz  force 


mdv 

dt 


5x6) 


(32) 
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The  magnetic  field  effects  can  be  divided  into  three  categories:  (1)  the  perpendicular  ©'ra¬ 
tion  about  the  magnetic  field  line  at  llj,  (2)  the  acceleration  along  the  field  by  the  mirror 
force 


Fa  =  -nTB, 


and  (3t  the  perpendicular  drift  across  the  magnetic  field,  given  by 


(33 1 


V.  =  cn 


Bx  VB 

eB- 


<341 


where  n  =  <  l/2lmvj_”,'B  is  the  magnetic  dipole  moment  of  the  particle.  The  parallel  accel¬ 
eration  of  the  particle  is 


aii  =  -(JKV® 


<3f» 


Note  that  aa  is  independent  of  the  mass  of  the  particle  and  of  the  magnetic  field  strength. 
It  depends  only  on  the  velocity  of  the  particle  and  the  scale  length  of  the  field.  For  parti¬ 
cles  ejected  perpendicular  to  the  field  at  v  =  3 km/sec  and  using 


V<nB  =  J 


(36) 


3 

Re 


(where  Rc  is  one  Earth  radius)  for  the  geomagnetic  dipole  field,  the  parallel  acceleration  is 


aj|  —  2m/sec^. 

The  magnitude  of  the  perpendicular  drift  velocity  is 


(37) 


=  h’. 

\  21 1, 


Vj^nB 


(38) 


which  is,  for  A  =  16  and  Z  =  1,  with  vj_  =  3  km/sec 


The  magnitude  of  Vj_  is  so  small  that  it  can  be  ignored  altogether. 

The  above  result  for  mirror  acceleration  holds  for  any  charged  particle,  including 
electrons,  ejected  at  3  km/sec.  The  thermal  speed  of  electrons, 


=  300  km  sec 

is  a  more  representative  velocity  of  the  electrons  in  the  exhaust.  The  parallel  acceleration 
for  this  velocity  is 

an  —  20km/'sec^  (41) 

while  the  perpendicular  drift  velocity 

Yl  -  0.3 cm/sec  (42) 

is  still  negligible.  For  Eq.  (41)  to  have  any  significance,  it  is  necessary  to  assume  that  the 
electrons  are  free  to  move  away  from  the  ions  in  the  exhaust. 

Plasma  processes  such  as  Alfven  critical  ionization  and  beam  plasma  discharge  are 
associated  with  unstable  growth  of  electrostatic  waves  which  accelerate  electrons  to  ener¬ 
gies  on  the  order  of  lOOeV.  In  general,  the  parallel  acceleration  of  electrons  is  given  by 

a||  =  1.8x  10^r(eV)r||fnBkm/sec^  (43) 

=  83t(eV)km/sec^ 

where  r  is  the  electron  energy  in  electron-volts.  Energetic  electrons  of  lOOeV  are  acceler¬ 
ated  at 

ajj  “  104km/sec2  (44) 

assuming  that  charge  separation  effects  are  not  important. 

The  rate  of  loss  of  energetic  electrons  by  parallel  mirror  acceleration  is  limited  by 
ambipolar  diffusion.  The  ambipolar  diffusion  coefficient  is 


(45) 


D  =  k(Te  +  T,) 
ntr(mewe  +  m1w,) 

where  n  is  the  neutral  number  density  (Chen,4  p.  160).  For  energetic  electrons  with 


The  time  scale  for  diffusion  losses  is 


where  L  is  a  typical  scale  length  in  the  problem.  Diffusion  is  important  only  if  tj  is  short 
compared  to  a  convection  time: 


The  ratio  of  diffusion  time  to  convection  time  is 

Td  LVmr 
tc  ~  w;  • 

Taking  L  —  r  and  using  Eq.  (5i,  we  have 


(471 


(48) 


Td  _  l.ixio3 

TC  r\  t  (eV) 


(49) 


The  distance  outside  of  which  diffusion  dominates  convection  (that  is,  T(J/tc  <  1  is) 


r  >  11  x  103  m 
r  2  -  *  m. 

\  c(eV) 


(50) 


So,  for  radii  greater  than  about  1  m  to  perhaps  100  m,  ambipolar  diffusion  of  energetic 
electrons  is  dominant  and  geomagnetic  mirror  acceleration  is  an  effective  process  for  re¬ 
moval  of  electrons  and  ions. 

The  effective  velocity  associated  with  removal  of  charged  particles  by  ambipolar  diffu¬ 
sion  can  be  estimated  by  equating  an  effective  mass  flux  to  the  diffusion  flux: 
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,.vCfr  =  dv,) 


r 

VI 


Near  the  spacecraft,  the  effective  diffusion  velocity  is  of  the  same  order  of  magnitude  as 
the  flow  velocity. 

A  charged  particle  injected  into  the  geomagnetic  field  has  a  bounce  motion  associated 
with  the  mirror  acceleration  because,  aside  from  slow  drifts,  the  particle  moves  along  a 
given  field  line  and  mirrors  at  each  end.  For  a  sufficiently  low  bounce  period,  particles 
mirroring  in  the  opposite  hemisphere  could  possibly  return  to  the  spacecraft,  affecting  the 
nature  of  the  interaction.  While  such  a  result  appears  to  be  unlikely  a  priori,  we  analyze  it 
in  the  following  manner:  The  bounce  period  is  roughly  the  local  radius  divided  by  the 
particle  velocity  (Schulz  and  Lanzerotti,0  p.  18): 


where  the  factor  .1.5  is  appropriate  for  a  mirror  latitude  of  23°,  typical  of  the  Shuttle  in 
low  inclination  orbit.  The  bounce  period  for  a  3km/sec  particle  is  about  7000 sec.  The 
bounce  time  of  a  nonrelativistic  electron  is 


I  tv 


7ii  ~  sec. 

11  \  rteV ) 


The  bounce  period  of  a  100 keV  electron  is  therefore  about  4  sec.  These  bounce  periods  are 
so  long  that  we  conclude  that  they  bear  no  likely  relation  to  the  Shuttle  interaction. 

An  important  question  in  our  analysis  is  the  depth  into  the  ionosphere  to  which  a 
charged  particle  penetrates  as  a  result  of  the  mirror  acceleration  for  a  given  velocity.  The 
reason  for  raising  this  question  is  to  investigate  possible  asymmetries  due  to  the  interac¬ 
tion  of  an  ionized  cloud  with  the  geomagnetic  field.  This  can  be  answered  approximately 
based  on  simple  one-dimensional  motion  arguments.  If  the  initial  velocity  and  constant 
acceleration  are  v|g  and  ajj,  then  the  penetration  depth  in  the  direction  along  the  field 
opposing  the  acceleration  is 


'el 


5.  Schulz,  M.,  and  Lanzerotti,  L.J.  (1974)  Particle  Diffusion  in  the  Radiation  Belts. 
Springer- Vcrlag,  Berlin. 
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Using  Eq.  (35)  for  the  parallel  acceleration,  the  penetration  depth  is  approximately 


vio^V  fnB 

_  ctn2« 

VfnB 


(56) 


where  a  is  the  initial  pitch  angle  (the  angle  between  the  initial  velocity  and  the  magnetic 
field).  It  is  clear  that  particles  with  any  significant  component  of  velocity  in  the  mirror 
direction  (toward  the  nadir)  will  travel  a  reasonable  fraction  of  an  earth  radius  in  that 
direction.  For  example,  particles  that  deviate  from  a  pitch  angle  of  90°  by  more  than  0.05° 
will  penetrate  more  than  1  km  toward  the  nadir. 

The  effect  of  large-scale  electric  fields  can  be  divided  into  two  categories:  (1)  parallel 
electric  fields,  E,|  =  B(B  E)/B^,  and  (2)  perpendicular  electric  fields,  Ej  =  E  -  E||.  The 
effect  of  a  parallel  electric  field  is  to  accelerate  the  particle  along  the  magnetic  field  at 
the  rate 


ZeE 

Am,, 


For  ions 


(57) 


aj  =»  9.6x  104^jE|j(V/m)km/scc2 


(58) 


and  for  electrons 


a|  =  1.8x  10®E||(V/m)km/sec^.  (59) 

For  purposes  of  illustration,  an  ion  with  A  =  16  and  Z  =  1  is  accelerated  in  a  field  Ej|  =  1 
mV/m  to  a  velocity  of 

Vjj  =  6km/scc  (60) 

in  lsec.  During  this  time  interval,  the  path  length,  starting  from  v||  =  0,  is 

s  =  3km.  (61 ) 
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The  corresponding  figures  for  electrons  are 


vj|  =  1.8  x  IQ*5  km/sec 


(62) 


and 


s  =  8.8  x  104km. 


(63) 


The  perpendicular  electric  field  causes  an  E  x  B  drift  which  is  the  same  for  electrons 
and  ions: 

V,  =  (64) 

x  B- 

This  drift  is  perpendicular  to  both  E  and  B;  its  magnitude  is 

Vj_  -  25Ei(V/m)  km/sec.  (65) 


The  velocity  in  a  1  mV/m  electric  field  is  0.25  km/sec. 

An  illustrative  parameter  in  determining  the  relative  effect  of  electric  and  magnetic 
accelerations  is  the  parallel  electric  field,  E||',  which  is  equivalent  to  the  mirror  accelera¬ 
tion.  This  electric  field  is  determined  by  the  equality 


^  -  (IF'V"8 

or,  for  ions, 

9TfnB 

E||'  =  Am„v^-7jj2e- 

=»  2.2  x  10"  8  |  V/m 


(66) 


(67) 


for  v  =  3km/scc.  For  electrons 

E,-  -  (t)r,(nB 


(68) 


17 


=  4.7x10  ^f(eV)  V/m 


Parallel  electric  fields  are  not  measured  to  be  greater  than  lmV/m  below  1500  km. u 
They  cannot  be  measured  at  all  below  about  lmV/m.  For  our  purposes,  we  would  like  to 
estimate  the  order  of  magnitude  of  the  maximum  electric  field  at  low  altitudes. 

In  a  static  atmosphere,  the  electric  field  is  determined  by 


VPj  =  njeE  +  njmig 

(69) 

VPe  =  -  neeE  +  nemeg 

(70) 

and 


TE  =  4rre(ni-ne) 


(71) 


where  g  is  gravitational  acceleration.  From  the  difference  of  Eqs.  (69)  and  (70) 


E  =  ne)[V(Pi~Pe)  -  nm,g| 


(72) 


we  have 


nkT  _  mjg 
(2neL)  P2e) 


(73) 


where  T  =  1000°  K  is  the  temperature  of  the  ionosphere.  It  is  apparent  that  E  cannot  be 
much  larger  than  the  largest  of  kT/(2eL)  and  mjg/(2c)  or  about  10“  ^  V/m.  Extremely 
small  electric  fields,  well  below  measurable  thresholds,  can  have  a  larger  effect  on  parallel 
acceleration  than  the  geomagnetic  mirror  force. 

This  pressure  balance  argument  ignores  motion.  The  electric  field  due  to  motion  of 
the  ionosphere  (due  to  frictional  coupling  of  the  ions  to  neutral  winds)  is 

E  =  -(J)vxB  (74) 


where  V  is  the  velocity  of  the  plasma  in  the  frame  of  the  earth.  Taking  a  velocity  of 


6.  Maynard,  N.  Private  communication. 
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0.1  km  sec  the  maximum  electric  field  is  of  the  order  of  0.002  V/m.  This  is  a  perpendicular 
electric  field  and  contributes  only  to  perpendicular  drifts.  This  electric  field  can  be  ne¬ 
glected  in  comparison  with  the  electric  field  induced  by  the  Shuttle  motion  (see  below). 

As  well  as  large-scale  electric  fields  intrinsic  to  the  plasma,  the  electric  field  due  to  the 
motion  of  the  spacecraft  should  be  mentioned.  If  the  electric  field  in  the  frame  of  the 
plasma  is  zero,  then  the  motion  of  the  plasma  as  measured  in  the  spacecraft  frame 


U  =  -Vs 


where  Vsc  is  the  spacecraft  velocity  measured  in  the  plasma  frame,  will  cause  an  electric 
field 


E  *  -(J)0xB. 


This  electric  field  is  just  such  as  to  cause  the  plasma  to  E  x  B  drift  at  the  velocity 


y  _  cExB 
B2 


in  the  spacecraft  frame.  The  Shuttle  velocity  is  U  =  7.3  km/sec  and  the  associated  electric 
field  is  0.2  V/m.  Any  ambient  perpendicular  electric  fields  would  have  to  be  of  order 
0.2  V/m  to  have  a  significant  effect. 

From  the  above  estimates,  we  learn  the  following: 

(1)  The  geomagnetic  mirror  force  is  insufficient  to  influence  the  behavior  of  ions  or 
electrons;  it  is  far  too  small.  In  addition,  we  may  make  the  following  observation:  any 
thermal  process  (critical  ionization,  for  example)  which  produced  large  numbers  of  ions 
would  be  expected  to  produce  an  isotropic  plasma,  at  least  as  far  as  the  direction  of  the 
field  is  concerned.  Consequently,  as  many  ions  or  electrons  would  travel  down  the  field 
line  as  up  because  the  mirror  force  cannot  turn  them  within  a  significant  fraction  of  an 
Earth  radius. 

<2i  The  behavior  of  ions  and  electrons  can  be  influenced  significantly  by  relatively  low 
level  (=0.1  mV/m)  parallel  electric  fields.  The  accelerations  and  resulting  velocities  are  on 
the  same  order  of  magnitude  as  velocities  resulting  from  the  motion  of  Shuttle. 

)3>  The  order  of  magnitude  of  the  parallel  electric  fields  required  is  much  smaller  than 
can  be  measured  experimentally  but  probably  larger  than  can  be  reasonably  expected  at 
low  altitudes. 

(4)  Ambient  perpendicular  electric  fields  are  not  expected  to  have  any  significant 
influence  on  the  evolution  of  ionized  contaminants  as  observed  in  the  frame  of  reference 
of  Shuttle. 

(5)  The  geomagnetic  bounce  time  is  probably  too  long  to  have  any  effect  on  short-lived 
( 1  sec)  Shuttle  exhaust  phenomena. 


•J’iWVvV 
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I.  MACROSCOPIC  PLASMA  BEHAVIOR 

The  interaction  of  the  plasma  and  electric  and  magnetic  fields  depends  on  the  density 
of  the  plasma.  The  behavior  of  tenuous  and  dense  plasmas  is  essentially  different.  Tenu¬ 
ous  plasmas,  for  which  the  rest  mass  energy  density  is  much  less  than  the  magnetic  en¬ 
ergy  density,  tend  to  behave  like  a  collection  of  single  particles  in  the  external  magnetic 
field.  On  the  contrary,  dense  plasmas,  for  which  the  rest  mass  energy  density  is  much 
greater  than  the  magnetic  field  energy  density,  tend  to  shield  any  external  electric  field 
and  move  across  magnetic  field  lines.  In  addition,  plasmas  whose  thermal  energy  density  i 
much  less  than  the  magnetic  energy  density  do  not  significantly  change  the  magnetic  field 
while  the  diamagnetic  effects  of  a  plasma  whose  thermal  energy  density  greatly  exceeds 
the  magnetic  energy  density  strongly  perturb  the  external  field. 

The  single  particle  behavior  is  governed  by  the  plasma  dielectric  constant 


Arrjic^ 

B2 


(see  Schmidt7).  This  is  the  ratio  of  the  rest  mass  energy  density  to  the  magnetic  energy 
density.  It  can  also  be  thought  of  as  c where  \\  is  the  (nonrelativistic)  Alfven 
speed,  or  as  the  square  of  the  ratio  of  the  plasma  frequency  to  the  ion  gyrofrequency: 


K  fl,2  ' 


Single  particle  behavior  is  determined  by  k  <  <  1.  If  k  >  >  1,  the  plasma  cloud  polarizes 
in  such  as  way  that  the  electric  field  in  the  center  of  mass  frame  (the  proper  frame  of  the 
plasma)  vanishes  and  the  momentum  of  the  cloud  is  conserved  independently  of  the  exter¬ 
nal  magnetic  field.  This  means  that  the  initial  momentum  of  the  plasma  cloud,  which  is 
determined  by  the  momentum  of  the  exhaust,  is  conserved. 

Using  the  values  of  the  exhaust  density  [Eq.  (4)1,  the  dielectric  constant  is  given  by 


k  =  3.2  x  10-°  f, 


where  f(r>  is  the  fractional  ionization  of  the  exhaust  and  r  is  in  cm.  The  distance  at  which 
k  —  1  is 


r  =  1.8  x  10''\  f  km 


7.  Schmidt,  George.  (1979)  Physics  of  High  Temperature  Plasma,  Academic  Press,  New 
York. 


VLN-xv<vw<yu 


If  there  is  strong  ionization  (for  example,  f  -  1),  then  the  plasma  polarizes  strongly  to  a 
distance  on  the  order  of  10°  km  and  does  not  follow  the  magnetic  field.  On  the  other  hand, 
single  particle  behavior,  for  which  the  particles  follow  the  ambient  magnetic  field  lines, 
past  about  10  m  requires 


f  -  10 
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Any  observation  which  showed  ions  travelling  along  the  ambient  magnetic  field  would 
suggest  that  the  ionization  state  was  extremely  weak.  The  actual  ion  density  implied  by 
F.q.  (82 1  is  about  lent  ~  'T 

The  diamagnetic  effect  of  the  plasma  on  the  magnetic  field  is  determined  by  the 
plasma  11: 


[1  = 


8-P 


(82) 


If  |1  <  <  1,  then  the  magnetic  field  is  essentially  given  by  the  ambient  field;  while  if 

|i  >  >  1,  then  the  field  is  strongly  perturbed,  essentially  killing  the  field  inside  the  plasma. 

Then  we  have 


11  -  4.5  x  1010  f  -- 
(Mr)J 


(84) 


where  M  is  the  Mach  number.  The  distance  for  which  (1  =  1  is 


r  = 


(85) 


If  the  plasma  is  highly  ionized,  the  magnetic  field  is  zero  to  perhaps  100  m  or  more,  while, 
if  the  fractional  ionization  is  on  the  order  of  10  ~  the  ambient  field  is  essentially  undis¬ 
turbed. 

These  considerations  lead  to  the  following  conclusions: 

1 1 1  A  sufficiently  fiigli  density  plasma  will  follow  the  initial  trajectory  of  the  exhaust 
for  a  large  distance.  If  the  exhaust  were  fully  ionized,  it  would  travel  a  perhaps  10° km 
before  it  was  tenuous  enough  to  follow  the  local  magnetic  field  lines. 

1 2)  A  sufficiently  low  density  plasma  would  follow  the  magnetic  field  immediately. 

This  leads  to  the  following  behavior  for  low  density  plasmas:  for  the  first  100  m,  they  are 
driven  radially  from  the  Shuttle  by  collisions  with  the  neutral  exhaust  (see  Section  1);  at 
about  100m,  they  follow  the  geomagnetic  field  and  consequently  move  away  from  the 
Shuttle  at  7 km  sec. 


(3)  A  fully  ionized  exhaust  kills  the  ambient  magnetic  held  until  about  10km  from 
Shuttle,  while  a  very  weakly  ionized  plasma  has  practically  no  effects. 


5.  NEUTRALIZATION  OF  POLARIZED  PLASMA  CLOUDS 

In  the  preceding  section,  it  was  implicitly  assumed  that  the  expanding,  highly  ionized 
cloud  remained  polarized  during  its  evolution.  In  fact,  the  polarization  state  of  the  cloud  is 
dependent  upon  the  surrounding  medium.  In  the  present  case,  ambient  electrons  are 
capable  in  principle  of  neutralizing  the  polarization  charge  density  which  shields  the  elec¬ 
tric  field  in  the  interior  of  the  cloud.  If  the  polarization  can  be  completely  neutralized,  the 
high  density  cloud  will  act  in  the  same  manner  as  a  low  density  cloud:  it  will  follow  the 
geomagnetic  field  lines.  The  purpose  of  this  section  is  to  inquire  about  the  time  required 
for  ambient  electrons  to  neutralize  plasma  clouds. 

A  rough  estimate  of  the  minimum  time  required  to  neutralize  a  plasma  cloud  is  the 
time  required  to  fill  the  cloud  with  ambient  thermal  electrons.  The  cloud  certainly  cannot 
be  neutralized  in  much  less  time.  The  time  rate  of  change,  due  to  thermal  electron  flux,  of 
the  charge  density  of  a  cloud  is 

dQ 

i,  =  -nwa  (86) 


where  w  is  the  electron  thermal  speed  and  a  is  the  cross  sectional  area  of  the  cloud  along 
the  field.  The  reason  for  taking  only  the  area  along  the  field  is  that  the  effective  electron 
velocity  across  the  field  is  the  Shuttle  velocity,  about  7  km/sec.  The  charge  in  the  cloud  is 


Q  =  fN 


(87) 


where  f  is  the  fractional  ionization  and  N  is  the  number  of  molecules  in  the  cloud.  The 
characteristic  time  for  neutralization  is  then 


r  =  -  Q/(dQ/dt) 


(88) 


=  fN 
nwa ' 


The  neutralization  time  is  seen  to  increase  with  the  ionization  and  depend  inversely  on 
the  flux  of  thermal  electrons.  Taking  N  =  5x  10”^,  n  —  lO^cm-**,  w  —  200km  sec,  and 
a  =  r^,  the  neutralization  time  is  seen  to  be  on  the  order  of 


t  —  2  x  10^2  f '  (cm )  sec. 


(89) 
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A  fully  ionized  cloud  cannot  be  neutralized  faster  than  about  10'*/r2sec.  For  example,  if 
the  extent  of  the  cloud  is  r  =  1km,  the  neutralization  time  is  lOOscc  or  more.  Neutraliza¬ 
tion  within  about  0.01  sec,  the  time  it  takes  the  exhaust  to  move  about  30m  at  3km/sec 
requires  f  to  be  less  than  about  10  "5.  We  conclude  that  weakly  ionized  clouds,  with  frac¬ 
tional  ionization  below  perhaps  10"  5  (and  likely  one  or  two  orders  of  magnitude  less 
because  of  the  qualitative  nature  of  the  argument),  can  be  neutralized  in  a  satisfactorily 
short  time.  Denser  clouds  continue  to  travel  across  the  field,  while  clouds  with  lower 
ionization  states  move  with  the  geomagnetic  field  almost  immediately. 

6.  DISCISSION 

From  the  discussion  in  each  section,  it  is  possible  to  develop  a  fairly  systematic  view 
on  the  expected  large-scale  evolution  of  Shuttle  exhaust.  Before  discussing  this  picture, 
however,  we  should  mention  that  a  number  of  interesting  and  important  topics  have  been 
omitted  entirely.  Among  these  are  the  source  of  the  ionization,  the  way  that  ionization 
develops  as  the  cloud  expands,  processes  such  as  two-stream  instabilities  associated  with 
streaming  of  ambient  electrons  through  the  ionized  cloud,  and  the  possibility  of  substan¬ 
tial  charge  separation  and  the  related  electric  fields  associated  with  the  collisions  of  ions 
and  electrons  with  the  neutral  exhaust. 

The  detailed  study  of  sources  of  ionization  is  crucial.  The  possibilities  range  from 
intrinsic  ionization  of  the  thermal  exhaust  to  Alfven  critical  velocity  ionization.  These  may 
be  expected  to  lead  to  low-level  and  high-level  ionization,  respectively.  From  what  has  been 
said,  it  may  not  be  possible  to  determine  low  levels  of  ionization  based  on  known  proper¬ 
ties  of  the  exhaust  because  of  the  extreme  sensitivity  of  the  ionization  level  to  the  temper¬ 
ature  of  the  exhaust  and  to  the  presence  of  contaminants  with  low  ionization  potentials. 
This  observation  underscores  the  importance  of  makingdirect  measurements.  In  fact, 
direct  measurements  of  the  radiation  intensity  and  perhaps  the  spectrum  of  ionized  ex¬ 
haust  are  crucial  to  determine  whether  the  ionization  level  is  low  or  high.  Without  such 
measurements,  and  an  accompanying  analysis  of  radiation  processes,  there  can  be  no 
certainty  of  even  the  approximate  nature  of  the  evolution  of  the  cloud.  The  measurements 
are  needed  to  guide  future  investigations,  essentially  to  establish  their  limits.  As  an  exam¬ 
ple,  suppose  that  measurements  establish  that  high-level  ionization  exists  in  Shuttle  ex¬ 
haust.  Then  the  nature  of  the  investigation  would  focus  on  the  Schmidt7  criterion  and 
neutralization  of  plasma  clouds.  On  the  other  hand,  measurements  establishing  low-level 
ionization  would  turn  the  investigation  toward  establishing  the  temperature  of  the  ex¬ 
haust  and  levels  of  contaminants  with  low  ionization  potentials. 

With  these  caveats  in  mind,  let  us  discuss  the  picture  of  the  exhaust  as  it  has  been 
developed  in  the  text.  The  driver  of  all  exhaust  phenomena  is  the  momentum  and  energy 
of  the  neutral  exhaust.  There  are  two  major  possibilities  for  ionization:  high  level  (with 
critical  ionization  as  a  likely  source)  and  low  level  (with  thermal  ionization  as  the  likely 
soured.  Up  to  a  distance  on  the  order  of  50  to  100  m,  the  evolution  of  the  exhaust  is  likely 
to  be  independent  of  the  ionization  level:  the  ions  are  dominated  by  collisions  with  the 
thermal  exhaust.  However,  near  100m,  the  exhaust  becomes  collisionless.  Past  this  transi¬ 
tion  region,  neither  ions  in  the  exhaust  nor  ambient  particles  are  stronglv  affected  by 
collisions.  This  has  two  major  implications: 


(1)  Since  neutrals  cannot  penetrate  deep  inside  the  exhaust  cone  inside  about  100  m, 
any  collisional  phenomena  such  as  chemical  reactions,  excitation,  or  ionization  by  ambi¬ 
ents  should  occur  outside  100  m.  Because  the  density  of  the  (neutral)  exhaust  is  larger 
near  100  m  than  at  greater  distances,  the  intensity  of  such  phenomena  will  probably  be 
greatest  near  this  distance. 

(2)  The  behavior  of  ions  inside  the  transition  distance  will  be  independent  of  ioniza¬ 
tion  level  inside  the  transition  distance  but  will  depend  heavily  on  the  ionization  level 
outside  it.  Low  density  ions  should  show  a  sharp  turn  at  the  transition,  where  they  are 
swept  away  at  up  to  7  km/sec  because  their  motion  is  dominated  by  the  geomagnetic  field. 
High  density  ions,  on  the  other  hand,  should  expand  radially  in  the  Shuttle  frame  of 
reference  to  a  large  distance  (many  kilometers  or  more)  before  they  can  be  affected  by  the 
field.  Between  these  two  limits,  there  should  be  a  continuum  of  trajectories.  For  moderate 
level  ionization,  above  a  fractional  ionization  of  10" 5  or  so,  the  neutralization  of  the  polar¬ 
ization  charge  of  the  plasma  cloud  cannot  be  accomplished  rapidly.  The  consequence  of 
partial  neutralization  is  that  the  cloud  should  be  partially  attached  to  the  magnetic  field 
and  be  swept  away  at  velocities  less  than  7  km/sec.  An  interesting  question  for  future 
investigation  is  the  nature  of  the  trajectory:  Is  there  still  a  sharp  break  near  100  m,  or  is 
there  a  smooth  trajectory?  The  answer  presumably  depends  on  the  details  of  the  neutral¬ 
ization  process,  which  have  not  been  treated  here.  The  qualitative  nature  of  the  paths  of 
low  and  high  density  ions  is  shown  in  Figure  3. 

The  behavior  of  a  highly  ionized  cloud  is  essentially  independent  of  ambient  electric 
and  magnetic  fields  for  a  large  distance.  As  such,  it  should  not,  for  example,  show  signs  of 
alignment  with  the  magnetic  field.  On  the  other  hand,  a  cloud  of  low  density  ions  can  he 
affected  by  these  fields.  It  can  be  accelerated,  for  example,  along  the  geomagnetic  field  by 
relatively  low-level  parallel  electric  fields  (  =  0.1  mV/m);  it  is  not  known  whether  such  large 
parallel  electric  fields  actually  exist  at  low  altitudes.  The  geomagnetic  mirror  force  is  not 
strong  enough  to  accelerate  ions  along  the  field.  As  a  result,  ion  clouds  cannot  be  expected 
to  show  signs  of  field  alignment.  The  sole  signature  of  the  presence  of  the  ambient  fields  is 
likely  to  be  the  sharp  break  of  low  density  ion  clouds  due  to  sweeping  away  of  the  cloud 
into  the  Shuttle  wake. 
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